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ABSTRACT 


We report on a LUCI/Large Binocular Telescope near-infrared (NIR) spectrum of PSO J030947.49+271757.31 (hereafter 
PSO J0309+27), the highest redshift blazar known to date (z~6.1). From the CIVA1549 broad emission line we found that 
PSO J0309+27 is powered by a 1.45*}82x10’M. supermassive black hole (SMBH) with a bolometric luminosity of ~8x10*° erg s! 
and an Eddington ratio equal to 0.444078. We also obtained new photometric observations with the Telescopio Nazionale Galileo in 
J and K bands to better constrain the NIR Spectral Energy Distribution of the source. Thanks to these observations, we were able to 
model the accretion disk and to derive an independent estimate of the black hole mass of PSO J0309+27, confirming the value inferred 
from the virial technique. The existence of such a massive SMBH just ~900 million years after the Big Bang challenges models of the 
earliest SMBH growth, especially if jetted Active Galactic Nuclei are associated to a highly spinning black hole as currently thought. 
Indeed, in a Eddington-limited accretion scenario and assuming a radiative efficiency of 0.3, typical of a fast rotating SMBH, a seed 
black hole of more than 10° Mo at z = 30 is required to reproduce the mass of PSO J0309+27 at redshift 6. This requirement suggests 
either earlier periods of rapid black hole growth with super-Eddington accretion and/or that only part of the released gravitational 


energy goes to heat the accretion disk and feed the black hole. 
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quasars: individual: PSO J030947.494+271757.31 


1. Introduction 


High redshift (z>6) Active Galactic Nuclei (AGNs) are direct 
probes of the Universe less than 1 Gyr after the Big Bang. 
These earliest AGNs are fundamental to study the early growth 
° of supermassive black holes (SMBHs, e.g. see Inayoshi et al. 
- 2020, for a recent review). An accurate determination of SMBH 
masses of z>6 AGNs is a prerequisite to fully understand 
SMBHs physics, demographics (e.g. black holes mass function), 
- « and relations with their host galaxies. With the so-called single- 
; = epoch (SE) method (e.g., Vestergaard & Peterson 2006; Vester- 
< gaard & Osmer 2009; Shen et al. 2008, 2011, 2019; Trakhtenbrot 
i 
ae) 
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& Netzer 2012) it has been possible to estimate the mass of the 
SMBHs hosted by high-z AGNSs (e.g., Jiang et al. 2007; Kurk 
et al. 2007; Wu et al. 2015; Mazzucchelli et al. 2017; Kim et 
al. 2018; Shen et al. 2019; Onoue et al. 2019; Yang et al. 2020; 
Andika et al. 2020; Wang et al. 2021). These studies have shown 
that high-redshift AGNs discovered to date are typically pow- 
ered by SMBHs more massive than 10°-10? Mo, comparable to 
the most massive black holes at any redshift. These discoveries 
indicate a fast and efficient growth of black holes that challenges 
the currently accepted theoretical model of SMBHs formation 
(e.g. Volonteri 2010, Latif & Ferrara 2016). The most popular 
scenarios to explain the mass assembly of several million solar 
mass black holes in the early Universe include direct collapse 


of massive gas clouds (e.g., Haehnelt & Rees 1993; Begelman 
et al. 2006; Latif & Schleicher 2015), the collapse of Popula- 
tion III stars (e.g., Bond 1984; Alvarez et al. 2009; Valiante et 
al. 2016), the co-action of dynamical processes, gas collapse and 
star formation (e.g., Devecchi & Volonteri 2009), or intense gas 
accretion 1n a super-Eddington phase (e.g., Alexander & Natara- 
jan 2014; Madau et al. 2014; Lupi et al. 2016; Pezzulli et al. 
2016; Volonteri et al. 2016). 

Even more challenging is the discovery of high-z massive 
SMBHs hosted in radio-loud!| (RL or jetted) AGNs. RL AGNs 
are characterized by the presence of two collimated relativistic 
jets of plasma emitted from the central SMBH and extended up 
to afew Mpc (e.g., see Blandford et al. 2019 for a recent review). 
The presence of relativistic jets is usually associated to a highly 
spinning accreting black hole (e.g., Blandford & Znajek 1977; 
Tchekhovskoy et al. 2011), which is expected to have a large 
radiation efficiency (7 ~0.3; e.g., Thorne et al. 1974) and, there- 
fore, for a given luminosity, a longer growth time with respect 
to black holes hosted by radio-quiet (RQ or non-jetted) AGNSs. 
Hence, since there is not enough time to accrete large masses 
(M,,, ~ 10? Mo) in a standard Eddington-limited accretion sce- 


' Here we consider an AGN to be RL if it has a radio loudness (R) 
larger than 10, with R defined as the ratio between the 5 GHz and 4400A 


rest frame flux densities, R = oils (Kellermann et al. 1989). 
4400A 
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nario, super-Eddington accretion events can be invoked to ex- 
plain the existence of these high-z jetted SMBHs (e.g., Begel- 
man & Volonteri 2017; Yang et al. 2020). Therefore, identifying 
and characterizing high-z RL AGNs provide a unique opportu- 
nity to study the role of jets in the accretion of SMBHs in the 
primordial Universe (e.g., Volonteri et al. 2015). 

If a RL AGN has its relativistic jets oriented along the line of 
sight, we classify it as a blazar (e.g., Urry & Padovani 1995; 
Padovani et al. 2017). Since the jet emission is strongly boosted 
and not obscured along the jet direction, the observed luminosity 
of blazars is usually very high, making these sources well visible 
up to very high-z. Although blazars represent a small fraction of 
RL AGNs, they are fundamental to ensure a reliable and com- 
plete census of the global population of the jetted AGNs and, 
therefore, to trace the evolution of the SMBHs across the cosmic 
time (e.g., Ajello et al. 2009; Ghisellini et al. 2010b; Sbarrato et 
al. 2015; Caccianiga et al. 2019, Ighina et al. 2021b; Diana et 
al. 2021). Indeed, from the space density of blazars it is possible 
to infer the space density of all the RL AGNs that have similar 
intrinsic physical properties. If we define as blazar a source ob- 
served within an angle equal to 1/I, where I is the bulk Lorentz 
factor of the emitting plasma, we expect to find Nez aGgns = 
Nolazars X 207, (e.g., Volonteri et al. 2011). Therefore, the dis- 
covery of high-z blazars ensures the census, free from obscu- 
ration effects, of early SMBHs and provides strong and critical 
constraints on the accretion mode, the mass and the spin of the 
first seed black holes (e.g., Kellerman 2016). 

Recently, we discovered the most distant blazar to date, 
PSO J030947.49+271757.31 (hereafter PSO J0309+27 at z~6.1; 
Belladitta et al. 2020). We are carrying out a multi-wavelength 
study on this source, from the radio (Spingola et al. 2020) to the 
X-ray band (Moretti et al. 2021, Ighina et al. 2022), in order 
to characterize the properties of this very distant jetted AGN at 
all wavelengths. Here we present new observations in the near- 
infrared (NIR) band. They consist on a Large Binocular Tele- 
scope (LBT) spectroscopic observation carried out to detect the 
CIVA1549 (hereafter CIV) emission line useful for the compu- 
tation of the central black hole mass, and on photometric obser- 
vations in J and K’ bands obtained at the Telescopio Nazionale 
Galileo (TING) to better constrain the Spectral Energy Distribu- 
tion (SED) of the source. 

The paper is structured as follows: in Section [2] we present the 
LBT and TNG NIR observations of PSO JO309+27; in Sec- 
tion [3] we report the data analysis and the results of our obser- 
vations (i.e., CIV line characterization, NIR magnitudes, black 
hole mass estimations); in Section [4] we discuss our results on 
the black hole mass in term of seed black hole growth; the con- 
clusions are reported in Section|5] 

Throughout the paper we use a flat ACDM cosmology, with Hg 
= 70 kms"! Mpc™!, Q,, = 0.3 and Qa = 0.7. All errors are re- 
ported at lo, unless otherwise specified. 


2. Spectroscopic and photometric observations 
2.1. LBT/LUCI 


The high-z nature of PSO J0309+27 (z~6.1) prevents the de- 
tection in the optical wavelength range of the CIVA1549 and/or 
the MgIIA2798 broad emission lines typically used for the black 
hole mass estimation. The latter, at the redshift of the object, 
falls in an atmospheric absorption band, hence it is not easily de- 
tectable by ground based telescopes. Therefore, we proposed an 
LBT Utility Camera in the Infrared (LUCI, Seifert et al. 2003) 
follow-up in order to extend the wavelength range in the NIR 
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band to detect the CIV broad emission line. The observation was 
carried out in a Director’s Discretionary Time program (program 
ID: DDT_2019B_3; PI: S. Belladitta) on 2019 December 2 and 
consisted of 12 exposures of 15 minutes each, with a long-slit 
of 1.2’ width, in nodding mode in the sequence ABBA, with a 
total integration time of 3 hours. The medium seeing throughout 
the night was 1.1” and the mean air mass was 1.2. We used the 
G200-zJ configuration for both LUCI1 and LUCI2, in order to 
cover the wavelength range from 0.9 to 1.2 um, where the CIV 
was expected to be found. The data reduction was performed 
at the Italian LBT Spectroscopic Reduction Center. Each spec- 
tral image was independently dark subtracted and flat-field cor- 
rected. Sky subtraction was done on 2D extracted, wavelength 
calibrated spectra. Wavelength calibration was obtained by us- 
ing several sky lines, reaching a rms of 0.33A on LUCII and of 
0.25A on LUCI2. The LBT/LUCI spectrum of PSO J0309+27 
is shown in Fig.|1} together with the LBT/MODS spectrum al- 
ready reported in Belladitta et al. (2020). Both the spectrum and 
the photometric points have been corrected for Galactic extinc- 
tion, using the extinction law provided by Fitzpatrick (1999), 
with a Ry = 3.1. The CIV emission line is clearly detected, 
and it can be used to compute the mass of the central SMBH 
of PSO J0309+27. 


2.2. TNG/NICS 


PSO J0309+27 was observed with J and K’ filters (central 2 = 
1.27 um and 2.12 um respectively) at the TNG during the night 
of February 12, 2021 (program ID: A42DDT4, PI: S. Belladitta) 
with the large field camera mounted on the Near Infrared Cam- 
era Spectrometer (NICS) instrument (Baffa et al. 2001) under 
excellent seeing conditions (FWHM~0.75”, see Table [ip. The 
observations consisted of 50-positions dithered-mosaic with a 
DIT of 1x60s (J-band) and 320s (K’-band) with a total in- 
tegration time of 50 minutes per band?! Flux calibration has 
been ensured by a short (30s) observation of the field of the 
AS1#|(RA = 05:57:07.5, DEC = 00:01:11) photometric stan- 
dard stars (ARNICA catalog; Hunt et al. 1998), just after the 
object acquisition. Standard data reduction, such as flat fielding, 
sky-subraction, cross-talk effect, image alignment and stacking, 
was performed with the Speedy Near-IR data Automatic reduc- 
tion Pipeline (SNA) properly optimized for NICS data. J and 
K’ images of PSO J0309+27 are reported in Fig. |2} The source 
appears point-like in both bands. The AB zero-point (ZP) mag- 
nitudes are reported in Table }1| along with the total exposure 
time and the seeing in the final mosaic, computed measuring the 
FWHM of reliable point-like objects by using the Image Reduc- 
tion and Analysis Facility (IRAF} task imexamine. 


3. Results and discussion 
3.1. CIV line width and luminosity 


We characterized the line width with both the Full Width at Half 
Maximum (FWHM) and the line dispersion (Gjjn¢, as defined in 


* Actually the final integration time in J band is 2940 s since one ex- 
posure was lost. 

3 There are four photometric stars in this field: AS13-0, AS13-1, 
AS13-2, AS13-3. However the AS13-0 is saturated, hence the photo- 
metric ZP was computed by using the remaining three. 


4 http: //www.tng.iac.es/news/2002/09/10/snap/index. 
http://ast.noao. edu/data/software 
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Fig. 1: LBT MODS and LUCI observed spectra of PSO JO309+27. Together with the optical lines already marked in the MODS spectrum in 
Belladitta et al. (2020), here the CIV11549 line is marked. On the upper x-axis the rest frame wavelengths are shown. 
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Fig. 2: 0.5’x0.5’ J and K’ cutout images of PSO J0309+27 taken with 
the LF camera mounted on the NICS instrument. Its optical position is 
marked with a red circle of 2”” of diameter. The two images are oriented 
with north up and east to the left. 


Peterson et al. 2004). We computed both by fitting the line pro- 
file, following different steps. First of all we de-redshifted the 
LBT/LUCI spectrum using a redshift of 6.063+0.003, which is 
based on the position of peak of the CIV emission line (Zcry). 
Then we linearly” fitted the continuum near the CIV line in two 
specific intervals (1445-1465A and 1670-1690A, see Fig.|3) free 


© AGNs usually show a power-law spectrum, but in short wavelenghts 
intervals the linear fit is a good approximation. 


Table 1: Details on TNG observations in J and K’ bands and 
PSO J0309+27 measured magnitudes. 


Filter = Acentral ZP seeing mag 

jum AB arcsec AB 

(1) (2) (3) (4) (5) 
J 1.27 23.25 0.76 20.81 + 0.06 
kK’ 2.12 23.92 0.75 20.93 + 0.08 


Notes. Col (1): NICS filter; Col (2): filter central wavelength; Col (3): 
photometric ZP in AB system; Col(4): Seeing (FWHM); Col(5): J and 
K’ AB magnitudes of the object. The relations to convert from Vega to 


AB systems are: Jag = Jvega + 0.91 and Ki, = Kvyeea + 1.85. 


from spectral features and spikes due to the background. On the 
pseudo-continuum subtracted spectrum we fitted the CIV broad 
emission line. Since the existence of a strong narrow component 
(produced by the Narrow Line Region) of the CIV line is con- 
troversial and difficult to detect (e.g., Wills et al. 1993; Corbin & 
Boroson 1996; Vestergaard 2002; Shen & Liu 2012), we did not 
include this component in the line fit. Moreover we decided not 
to include the Fell features, because, as mentioned in previous 
studies (e.g., Shen et al. 2011; Trakhtenbrot & Netzer 2012; Zuo 
et al. 2020), the contribution from Fell around the CIV line is 
expected to be small. From Fig. [3] it is clear that two Gaussian 
functions are necessary to reproduce the broad emission line pro- 
file properly. Indeed, it has been already demonstrated in several 
works (e.g., Laor et al. 1994; Shen et al. 2008; Tang et al. 2012) 
that the CIV broad line is usually well described by a multiple 
Gaussian profile and not by a single Gaussian function. The two 
Gaussian components shown in Fig. [3] have the following char- 
acteristics: the larger Gaussian is centered at ~ 1541A and has 
a FWHM of ~ 9520 km s7!; the narrower component, at A = 
1549 A, has a FWHM of ~ 1330 km s7!. In Table|2|we reported 
the best fit parameters for the total CIV emission line (1.e., those 
derived from the sum of the two Gaussian components). Besides 
the FWHM (15.3 A) and the Cjine (19.7 A), from the best fit 
model we also measured the rest-frame equivalent width (REW) 
and the line luminosity (Lcry = 47D7 Fcrv). The uncertainties on 
these values were evaluated through a Monte Carlo method (e.g., 
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Shen et al. 2011; Shen & Liu 2012; Raiteri et al. 2020; Zuo et al. 
2020, Diana et al. 2021). Each wavelength of the best fit model 
was randomly perturbed for 1000 times, according to a Gaussian 
distribution of the mean rms of the spectra computed underneath 
the CIV line on the pseudo-continuum subtracted spectrum. In 
this way we obtained 1000 different mock spectra of the line pro- 
file, from which we measured the line properties with the same 
procedure used on the real data. We computed the distributions 
of FWHM, Ojine, REW and Fery for these 1000 simulated spec- 
tra, and the interval that contains 68% of the data in these dis- 
tributions was taken as the statistical uncertainty on the best fit 
values. 


flux rest [10~2” erg/s/cm?/A] 





1450 


1500 1550 


A rest [A] 


1600 1650 1700 


Fig. 3: Double Gaussian fit of the CIV emission line. The rest frame 
LBT/LUCI continuum-subtracted spectrum is reported in grey, the two 
Gaussian functions in red and the sum of the two in blue. The larger 
Gaussian (solid red line; 2 = 1541A) has a FWHM of ~9520 km s7!; 
instead the narrower component (dotted red line; 2 ~ 1549A) has a 
FWHM of ~1330 km s7!. The intervals for the continuum selection 
are indicated with dashed green vertical lines. 


3.2. CIV line blueshift and asymmetry 


The CIV emission line is known to show asymmetry and to be 
blueshifted with respect to low ionization lines (e.g., Gaskell 
1982; Richards et al. 2011; Coatman et al. 2017; Vietri et al. 
2018; Zuo et al. 2020), independently from the source orien- 
tation (e.g., Kimball et al. 2011; Runnoe et al. 2014). These 
characteristics suggest that the CIV clouds are affected by non 
gravitational effects, such as outflows, most likely originated in 
disk winds. Large CIV blueshifts indicate that non-virial mo- 
tions have a significant effect on the observed emission velocity 
profile. To date, the largest blueshifts (>3000 km s~') have been 
discovered in the so-called weak emission line quasars (WELQs, 
Diamond-Stanic et al. 2009), which exhibit a REW < 10A anda 
strongly asymmetric line profile (see e.g., Vietri et al. 2018 and 
reference therein). With a REW of ~25A PSO J0309+27 does 
not belong to this quasar population. Therefore, we will not ex- 
pect to find a high value of blueshift for our source. 

Moreover several works in the literature (e.g., Shen & Liu 2012; 
Coatman et al. 2016) found that for small values of FWHM and 
Cline (<5000 km s~!), a low value of blueshift (Av<2000 km s~!) 
is usually observed. 
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We computed the CIV line blueshift of PSO J0309+27 fol- 
lowing the equation of Coatman et al. (2017): Av (km s7!) = 


oe where c is the speed of light, 1549.48A is the rest 


frame wavelength for the CIV and Ajgi is the line centroid’ 
The value of the estimated Ajyqj¢ and Av is reported in Table |2 
The blueshift value is smaller than 2000 km s~!, as we expected, 
which is an indication that the outflows component is weak with 
respect to the emission of virialized gas. This allows us to infer 
that the virial black hole mass of PSO JO309+27 computed in 
Sect.[3.4.1]should not be strongly affected by blueshift effects. 


3.3. Object magnitude in J and K bands 


Once measured the ZP thanks to the photometric stars, we de- 
rived the J and K’ magnitude of PSO J0309+27 by using the 
IRAF aperture photometry package gphot and the extinction 
curve of the observational site. The aperture size was chosen 
to be three times as large as the seeing. J and K’ AB magni- 
tude of PSO J0309+27 are reported in Table [i] The conversion 
factor (in Vega system) to switch from K’ to K magnitude has 
been computed by convolving K and K’ NICS filter transmission 
curves with an AO stellar template. We found: K’ = K + 0.159. 
Therefore, we obtained: K (AB) = 20.77 + 0.08. The extracted 
J magnitude is consistent with that found for PSO JO309+27 in 
the UKIRT Hemisphere Survey (UHS, Dye et al. 2018): 19.51 + 
0.38 (Vega) = 20.42 + 0.38 (AB). 


3.4. Black hole mass estimation 


We computed the central black hole mass (Mpy) of 
PSO J0309+27 following two different and independent meth- 
ods. The first method is the commonly used virial approach (the 
single epoch, SE, method) and the second is based on the mod- 
eling of the accretion disk emission. 


3.4.1. Single epoch mass 


The SE approach is the most used and reliable method to com- 
pute black hole masses of Type I un-obscured AGNs. Although 
some works have questioned the reliability of CIV as a good 
virial mass indicator (e.g., Sulentic et al. 2007; Shen & Liu 2012; 
Trakhtenbrot & Netzer 2012) due to its observed blueward asym- 
metry and velocity shifts of the line profile, other authors have 
demonstrated that there are no large inconsistencies between the 
SE Mpgu computed from CIV and Balmer lines (e.g., Vestergaard 
& Peterson 2006; Greene et al. 2010; Assef et al. 2011, Dalla 
Bonta et al. 2020). 

To compute the black hole mass of PSO JO309+27 we followed 
the scaling relation of Vestergaard & Peterson (2006) based on 
the o lind’ 


o(km/s) 
10°km/s 


A LAr ssoa 


0.53 
—___—— 1 
10“erg/ - ) 


Mgu = 10°”? x ( ) x 


We used this relation because Denney et al. (2013) and Dalla 
Bonta et al. (2020) find better agreement between CIV-based 


’ The line centroid is defined as the wavelength that bisect the line in 

two equal part. We used the definition of Dalla Bonta et al. (2020): Anaiy 
o_o where P(A) is the line profile. 

§ Tn Appendix|A. 1}we reported the Mgy derived by using the FWHM, 

to better facilitate the comparison with SMBH masses of other high-z 

AGNs in the literature. 
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Table 2: Best fit parameters for the total CIV broad emission line. 


FWHM Cline REW Forv Lev Anaif Av 
km s7! km s7! A 10°'Sergs''cm-? 10“ erg s“! A km s7! 
(1) (2) (3) (4) (5) (6) (7) 
29607137 3R 1571 25,3707 8567024 35470101542, 16°09S [429+ 





Notes. Col (1) and (2): line width in term of FWHM and oj; Col (3): rest-frame equivalent width; Col (4) and (5): line flux and line luminosity; 


Col (6): line centroid; Col (7): line blueshift. 


and Hg-based mass estimates by using oj,- rather than FWHM, 
in particular when high quality spectra are used. Moreover we 
did not use black hole mass estimators that correct the effect of 
the CIV line blueshift (e.g., Coatman et al. 2017) for the follow- 
ing reasons: 1) the Ojjn¢ parameter is relatively insensitive to the 
blueshif(}| (e.g., Coatman et al. 2017, Dalla Bonta et al. 2020); 
2) Coatman et al. (2017) suggest to use these estimators when 
the blueshift value is larger than 3000 km s~! (the values mea- 
sured for PSO J0309+27 is not so high, see Sect. (3.2); 3) the 
application of the blueshift correction factor, calibrated on z<4 
AGNs may be inappropriate for sources at higher redshifts (e.g., 
Park et al. 2017; Mejia-Restrepo et al. 2018; Kim et al. 2018). 
The continuum luminosity at 1350A (ALA,359,) has been com- 
puted directly from the PSI y point (rest-frame wavelength = 
~1370A): ALj,,., = 2.49+0.32x 10" erg 1. 


Therefore, from Eq. ‘I]we obtained a Mgy equal to 1.457753 x 10° 
Mo. The reported uncertainty, derived by propagating the errors 
of the CIV line width and of the monochromatic continuum lu- 
minosity, is purely statistical. By taking into account the large 
intrinsic scatter of the CIV relation of the SE method (~0.36 dex, 
e.g., Vestergaard & Peterson 2006; Denney et al. 2012; Jun et al. 
2017), that dominates the overall Mgqy error, we obtained a black 
hole mass of 1.45*)8:x10°Mo. 

Since PSO J0309+27 is an object observed under a small view- 
ing angle, it is important to take into account the potential prob- 
lematics related to the use of the SE method on this type of 
source. First, since the Broad Line Region (BLR) may not be 
isotropic the resulting black hole mass could be systematically 
underestimated in objects observed face-on. Indeed, there is a 
large consensus in the literature (e.g., McLure & Dunlop 2002; 
Decarli et al. 2008, 2011) concerning to a disk-like structure of 
the BLR. However, it is not already clear if the width of the broad 
emission lines (including CIV) depends on the orientation. For 
instance, Runnoe et al. (2014, and reference therein) found this 
dependence in RL AGNs for H{, but not for the CIV line. Sim- 
ilarly, Fine et al. (2011) in a sample of RL AGNSs did not find a 
correlation between the line width of the CIV and the AGN ori- 
entation. These authors concluded that the high ionization lines 
are produced in isotropic inner parts of the BLR. 

A second potential issue connected with the orientation of the 
source is the fact that the AGNs continuum luminosity could 
be contaminated by the relativistic jet. Therefore, a continuum- 
luminosity based relationship may lead to a mass overestimate 
(e.g., Decarli et al. 2011). Moreover, we also have to take into 
account the possible anisotropy of the continuum emission pro- 
duced by the accretion disk. Therefore, the observed continuum 
luminosity is higher for a source viewed face-on (e.g., Calderone 
et al. 2013). This effect could lead to overestimate mass, since 
the SE relations are empirically calibrated on type-I AGNSs ran- 
domly oriented (with an expected mean angle of ~30°). How- 


° The blueshift correction for black hole mass estimators has been cal- 
ibrated only for the FWHM parameter (e.g., Coatman et al. 2017). 


ever, there are no evidence of the presence of these potential bias 
for oriented RL AGNs, as it has been recently demonstrated by 
Diana et al. (2021, MNRAS, submitted). In this work the authors 
did not find a significant difference between the ratio between the 
CIV line luminosity (which is not affected by the beaming) and 
the continuum luminosity at 1350A (which could be affected by 
the beaming) of a sample of ~380 blazars with that of the RQ 
AGNs of the sample of Shen et al. (2011), for which the beam- 
ing 1s not present and that are, on average, observed at different 
angle compared to blazars. We computed this luminosity ratio 
(R=Lj3594/Lcrv) also for PSO J0309+27, finding that it is at lo 
from the mean value of the sample of Shen et al. (2011). This 
allows us to infer that the peculiar orientation of PSO J0309+27 
does not affect the observed continuum emission and, hence, the 
derived SE black hole mass. 


3.4.2. Accretion disk model 


To verify the presence of any possible bias on the calculated SE 
masses we used an independent technique based on the accretion 
disk emission (e.g., Sbarrato et al. 2012; Calderone et al. 2013; 
Ghisellini et al. 2015; Belladitta et al. 2019; Paliya et al. 2020; 
Diana et al. 2021, MNRAS, submitted). This technique assumes 
that the optical/UV continuum emission of the AGN is produced 
by an optically thick, geometrically thin accretion disk (AD) that 
emits according to the Shakura & Sunyaev (1973, SS73) model. 
The SS73 assumes a non-spinning!®| black hole (i.e., the effi- 
ciency of the accretion process, 7, 1s ~0.1), surrounded by an 
AD divided in rings that emit as black bodies. The total disk lu- 
minosity is therefore a superposition of black body spectra with 
the form: 


Rout 
L(V, Mgy, M)dy = 4 { RB,[T(R, Mgy, M)|dvdR (2) 


Rin 


where M is the mass accretion rate, R is the distance from the 
central engine and B,[7T(R, Mgx, M)]dyv is the Planck’s spec- 
trum. The temperature of the single black body emission de- 
pends on the distance from the central black hole, with the high- 
est temperature emitted at the disk inner radius, which corre- 
sponds to the innermost stable circular orbit: Rin = 3Rschw = 
6G Mpy/c? (e.g., Misner et al. 2017 and reference therein). A 
large black hole mass implies a larger disk inner radius and 
hence a lower emitting temperature (as T « May, an Therefore 
the accretion disk radiation peaks at lower frequency compared 
to smaller masses. 

With these assumptions it is possible to derive the values of M,,, 


10 Usually RL AGNs are associated to spinning black holes. However 
the assumption of a non-spinning black hole is justified also in the 
case of PSO J0309+27 by the results of Campitiello et al. (2018) who 
found an equivalence between the accretion disk fit with a SS73 model 
and a KerrBB model with spin~0.8 observed face-on (as expected for 
blazars). 
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and of the accretion rate (M), that are free parameters of the 
SS73 model, simply by fitting the optical/UV data points that do 
not suffer from HI absorption, 1.e. those at frequency lower than 
the Ly-a@ line. With the main goal of testing the SE masses, we 
decided to apply the AD method to PSO J0309+27 even if the 
optical/UV spectral range is not well sampled (see Fig. |4p, since 
we have only three photometric points not affected by HI ab- 
sorption: the y data from PS1 and the data at J and K bands from 
the TNG follow-up. In addition, the continuum of the LBT spec- 
trum has been used as a guide line for the modeling. By using 
this method we found that PSO J0309+27 can be described by a 
black hole of Mgy = 6.97;310* Mo. To reduce the uncertainty 
interval on the estimated aes we have used the expected value 
of the peak luminosity (Lypeax) of the accretion disk as a further 
constraint. Lyeqx can be inferred from Leyy and/or from AL... . 
by using the relations found in Calderone et al. (2013, equations 
(2) and (5) of this paper) and by taking into account the inclina- 
tion of the source (i.e, the fact that PSO J0309+27 is a blazar, 
§ ~ 0°). This constraint allowed us to infer that the value of the 
black hole mass of the source cannot be lower than 5.2x10® Mo. 
Therefore the final black hole mass of PSO JO309+27 computed 
with the accretion disk method is: Mgy = 6.97)7x10® Mo (see 
Fig.|4). This value is consistent with that ébtained from the virial 
method, considering the total uncertainties. 


— MODS +LUCI spectra 
—— SS73: Log(Msxy/Mo) = 8.83 
@  PS1 and TNG photometric points 


log(vL,) [erg/s] 





14.8 14.9 15.0 15.1 15.2 15.3 15.4 15.5 15.6 15.7 
log(v) rest-frame [Hz] 


Fig. 4: Accretion disk model of PSO J0309+27. The SS73 model that 
best represents our data is shown in red. Orange dashed line and green 
dashed-dotted line represent the model with the minimun and maximun 
value of Mgy respectively (5.2x10® and 8.8x10° Mo). In black the opti- 
cal and NIR (smoothed) spectra are reported. The PS1 y and the TNG J 
and K photometric data are shown in light blue. In grey we highlighted 
the photometric point and the part of the spectrum affected by HI ab- 
sorption (therefore not used for the disk modeling). 


3.5. Bolometric Luminosity and Eddington ratio 


Using the value of the SE black hole mass, we derived the Ed- 
dington luminosity (Lggq) and the Eddington ratio (Aggqg). The 
latter quantifies how fast the accretion rate is with respect to the 
Eddington limit. To compute Aggqa we first estimated the bolo- 
metric luminosity (Lg,;) of PSO JO309+27 using a bolometric 
correction (e.g. Richards et al. 2006): Logo: = Ly3594 X Kooi. In 
this case we used the bolometric correction factor of Shen et al. 
(2008): K,o; = 3.8141.26. However we have to remind that this 
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Kgoi 18 calibrated empirically over RL and RQ Type I AGNs, 
with a mean expected angle of 30°. Since PSO J0309+27 is a 
source seen under a small viewing angle (@ ~ 0°) and that the 
continuum emission from the disk is not isotropic, as already 
mentioned in Sect. [3.4.1] we have to take into account the ex- 
pected inclination factor (i= os = = 1.15) to compute the intrin- 
sic bolometric luminosity. This led to a final estimate of: Lyo; 
= 8.22+3.70x10* erg s~!. We obtained a similar value for Ly, 
(~8x10* erg s~!) by using the non linear relation between Ly, 
and L359, of Runnoe et al. (2012). Then we computed Agag as 
as the ratio between the bolometric luminosity (1.e. that includes 
the optical/UV radiation of the accretion disk, the emission re- 
processed by the molecular torus 
tion) and the Eddington luminosity derived from both the Gjn¢ 
and the FWHM"| HY Tn the following sections we always use the 
Arag estimated from the o7;,-, which we considered the best BH 
mass estimator. The obtained value is: Agaa = 0.4473-{8. The un- 
certainty already takes into consideration both the statistical er- 
ror on the virial mass and the intrinsic scatter of the SE relation 
(~0.36 dex). 

The values of the SE Mpy and of Agga of PSO JO309+27 are in 
line with those derived for RL and RQ AGNs discovered at sim- 
ilar redshift (2=5.5-6.5]"| Broadly speaking, PSO J0309+27 is 
fully consistent with a typical z~6 AGN and shows no evidence 
of peculiarities associated with its relativistic beamed jet. How- 
ever the similarity of masses and Eddington ratio could be likely 
a consequence of a selection bias, as all these high-z sources 
have been selected from similar optical/IR surveys. 

If we consider, for the Aggg computation, only the luminosity of 
the accretion disk, which is Lgjsk ~ Lyoi/2 (e.g., Calderone et al. 
2013), we obtain an Eddington ratio equal to 0. 22404 a 








4. Implication on the early SMBH growth 


Accurate measurements of black hole masses and Eddington ra- 
tios of high-z AGNs help us constraining formation scenarios of 
the first seed black holes. Moreover, high-z RL AGNs provide 
a unique opportunity to study the role of jets in the accretion of 
early SMBHs (e.g., Volonteri et al. 2015). In particular, as al- 
ready mentioned, the discovery of jetted AGNs in the early Uni- 
verse (z > 5) represent a serious challenge to our understanding 
of black hole growth, especially if the presence of the jet is as- 
sociated with a rapidly spinning black hole, which is expected to 
have a large radiation efficiency and, therefore, a longer growth 
time with respect to black holes hosted by RQ AGNs. 
Assuming that the black hole seed grows at a constant Edding- 
ton ratio during the entire accretion process (e.g., Shapiro 2005; 
Volonteri & Rees 2005), the evolution of the Mpy with time is 
directly proportional to the mass itself, resulting in an exponen- 
tial growth from the initial mass (May seed): 


rowth 
Msu, seed — = Meu x exp(- oan (3) 
where tg;owrn 18 the time during which the black hole accretes and 
T 1s the e-folding timescale: 


] 
ea) mre ee | 
— 7) 


G 4 
ia ) Gyr (4) 


Fact(M, t) 


'l See Appendix[A] for details. 

'2 SE black hole masses of 5.5<z<6.5 RQ and RL AGNs reported in the 
literature varies from 8x10’ M, to 1x10!°M,; the value of Agua varies 
from 0.03 to 1.3; Jiang et al. 2007; Willott et al. 2010b; De Rosa et al. 
2011; Mazzucchelli et al. 2017; Eilers et al. 2018; Kim et al. 2018; Shen 
et al. 2019; Onoue et al. 2019; Andika et al. 2020. 
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where fact 18 the duty cycle of the black hole, i.e. the mass and 
redshift dependent fraction of time when the black hole 1s active 
(fact(M, t)=[0,1]). Therefore, from Eq. [3] we can derive the ini- 
tial mass of the black hole seed required to observe the mass of 
PSO J0309+27 at z~6 (0.922 Gyr after the Big Bang). Figure|5] 
shows the estimated growth history of PSO JO309+27 according 
to this model. These results depend on the assumptions made, 
i.e. the redshift of the seed formation, the accretion rate, the 
radiative efficiency, the value of fa.;. We traced the mass back 
to z=30, when the first stars and galaxies are thought to have 
formed (e.g., Bromm & Larson 2004; Bromm & Yoshida 2011). 
The value of the duty cycle (f,<;) is assumed to be equal to 1, 
meaning that the AGN has been active for all the time. Then we 
considered different values for 7 and Agag. We assumed that the 
seed black hole accretes constantly with the observed Eddington 
ratio (Agaq = 0.44) or with the value of 1, the maximum value al- 
lowed in an Eddington limited accretion scenario. The efficiency 
parameter, instead, is believed to depend on black hole spin (e.g., 
King & Pringle 2006) and can be as high as ~30-40% in case of 
spinning black hole (e.g., Thorne 1974; Reynolds 2014). Cur- 
rent semi-analytical models place only weak constraints on the 
spin values for AGNs at z>5, which depend on the gas accre- 
tion mode, the morphology of the host galaxy and black hole 
mass (e.g., Sesana et al. 2014). Therefore, as no stringent con- 
straints on black hole spin are reported to date for high-z AGNs, 
we assumed both an efficiency of 0.1 (typical of not rapidly spin- 
ning black holes) and 0.3 (typical of Kerr black holes). Figure|5| 


ke PSO~J0309+27 z=6.1 


n= 0.3 Agagg=1 


n= 0.1 Acad obs 


Direct collapse 


Stellar dynamical processes 


Pop-lll stars 





108 109 
Time since Big Bang (years) 


Fig. 5: Estimated growth history of PSO JO309+27 (red star). Solid 
lines represent the best fit cases, under which the corresponding Aga 
and 7 are reported. The shaded horizontal regions correspond to the 
expected mass ranges of Pop III remnants BHs (M,..g<107Mo, grey), 
stellar dynamical processes (M,..~210°-10*Mo, cyan) and direct col- 
lapse BHs (Myeeg~10*-10°Mo, dark cyan). The values of the different 
seed black holes are taken from Valiante et al. (2016). 


shows that only a scenario of 7=0.1 can reproduce a theoretically 
accepted seed mass. Scenarios of higher efficiency (7=0.3), in- 
stead, would require more massive seeds (M seeq> 10°M.) as pro- 
genitors of PSO J0309+27. These expected seeds are even more 
massive than what direct collapse models predict (Myeeq~10*- 
10°Mo, Latif & Ferrara 2016). 

This result suggests that such high values of efficiency are 
probably not realistic, not even for RL AGNs. Alternatively, 


super-Eddington accretion episodes must occur for a signifi- 
cant fraction of the growth time. To date, there have been no 
clear examples of such super-Eddington SMBHs at z>6, al- 
though this scenario has been suggested for J1205—O000, a 
mildly obscured AGN at z=6.699 (Onoue et al. 2019), for 
PSO J006+39 at z=6.621 (Tang et al. 2019), two RQ AGNs, 
and for PSO J172+18, the most distant RL AGN ever discov- 
ered (z=6.8, Bafados et al. 2021). Super-Eddington accretion 
episodes are often taken into consideration also for the growth 
of black holes hosted in the z>7 RQ AGNs discovered so far 
(e.g., Banados et al. 2018a; Wang et al. 2021). It has been sug- 
gested that maintaining super-Eddington accretion might be pos- 
sible in specific environments (e.g., dust-obscured AGNs with 
strong winds or gas rich AGNs; Kim et al. 2015; Kubota & Done 
2019; Moffat 2020), but whether or not this type of accretion is 
sustainable remains an important open question for the growth 
of both RL and RQ AGNs. 

Another possible solution taken into consideration for the growth 
of black holes in RL AGNSs has been proposed by Jolley & Kun- 
zic (2008), Jolley et al. (2009) and Ghisellini et al. (2010a). 
These authors proposed that when a jet is present, not all the 
gravitational energy of the infalling matter is transformed into 
heat and then radiation, but can be transformed in other forms, 
such as amplifying the magnetic field energy of the inner disk, a 
necessary ingredient for launching the jet (Blandford & Znajek 
1977). In this case the total efficiency of the accretion process 
can be equal to 0.3, but only a fraction of it (7g, 1.e. the radiation 
efficiency) goes to heat the disk, while the rest (1-77) amplifies 
the magnetic field necessary to launch the jet. Therefore, disk lu- 
minosity becomes Eddington limited for a larger accretion rate, 
making the black hole growing faster. 

Cosmological simulations of seed black holes growth (e.g., Di 
Matteo et al. 2008; Alexander & Hickox 2012; Feng et al. 2014a) 
are fundamental to understand what are the main ingredients of 
the black hole seed evolution besides the accretion process. In 
particular, mergers and AGN feedback could be taken into ac- 
count for understanding the evolution of SMBHs hosted in RL 
AGNs. Indeed, RL AGNs are commonly found in rich environ- 
ments at different cosmic epochs (from redshift 0.5 to z=5.8; 
e.g., Pentericci et al. 2000; Venemans et al. 2002, 2004; Zheng 
et al. 2006; Hatch et al. 2014). Theoretical models strongly sup- 
port a preferential over-dense environment around RL AGNs 
(e.g., Orsi et al. 2016; Izquierdo- Villalba et al. 2018) and sim- 
ilar conclusions have been also found by studying the RL AGNs 
level of clustering with cosmic times (e.g., Magliocchetti et al. 
2004; Retana-Montenegro & RO6ttgering 2017). All these re- 
sults suggest that the presence of a relativistic jet may indeed 
be preferentially triggered in dense environments (1.e., in proto- 
clusters), Where frequent mergers between star-forming galax- 
ies help to increase the mass and spin of the SMBHs (e.g., 
Hatch et al. 2014). In this context a study of the environment 
of PSO J0309+27 could be crucial to better understand its black 
hole growth. 


5. Summary and conclusion 


In this paper we have reported new photometric and spectro- 
scopic observations in the NIR band of of PSO J0309+27, the 
most distant blazar discovered thus far. From a LUCI/LBT spec- 
troscopic observation we have detected the CIVA1549 broad 
emission line, which allowed us to compute the mass of the 
SMBH hosted by the source. By parameterizing the CIV line 
width with the oj,- and by using the SE method we estimated 
a mass for the central SMBH of 1.45*582x10" Mo. Moreover, 
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thanks to a dedicated follow-up with the TNG in J and K’ bands 
we also better constrain the NIR SED of the source, allowing us 
to derive and independent estimate of the SMBH mass using a 
method based on the accretion disk emission. The agreement of 
the two results supports the reliability of our estimate. 

The value of the black hole mass, the Eddington ratio and the 
bolometric luminosity of PSO J0309+27 are in line with those 
of other RQ and RL AGNs at similar redshift. However, to fully 
understand if the high-z blazar population is different from the 
RQ and/or RL ones in term of black hole mass and Aggq, a larger 
and statistically complete sample of blazars at the highest red- 
shift is required. 

Finally, we have computed the mass of the seed black hole 
required to reproduce the mass of the SMBH hosted by 
PSO J0309+27, using a simple model for the SMBH growth. We 
found that to obtain a reasonable (i.e., predicted by the models) 
seed black hole the efficiency of the accretion process can not be 
as high as 0.3, as expected for a SMBH hosted by a RL AGN. 
A high efficiency of 0.3 could be possible if super-Eddington 
accretion episodes are taken into account during the black hole 
growth or only a part of the released gravitational energy of the 
infalling matter is used to heat the accretion disk. 

Future studies on the environment of PSO J0309+27 will be use- 
ful to better understand its growth and evolution. 
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Appendix A: PSO J0309+27 black hole mass and 
parameters computed from FWHM 


In this work we have used the SMBH mass based on the Oj. 
parameter as the best Mgy estimator. Here we report the black 
hole mass, and the parameter related to it, computed from the 
FWHM, for a direct comparison with the literature. 

The single epoch scaling relation of Vestergaard & Peterson 
(2006) based on the FWHM is: 


FWHM(km/s).,__, ALAy3504 
10°km/s 10¥erg/s 


The value of the FWHM of the CIV line has been computed 
in Sect. and it is reported in Tab. |2} The luminosity of the 
AGN continuum is reported in Sect. [3.4] FromJA.1]we derived a 
black hole mass of 7.AT*S-3x108Mo, that is consistent with that 
obtained by ojjne and with the AD method. Also on this estimate 
we have to take into account an intrinsic scatter of ~0.4 dex. 
Since the value of the blueshift of PSO JO309+27 is smaller than 
3000 km s7! (see Sect. (3.1, by using the equation of Coatman 
et al. (2017) that correct the virial mass for the blueshift effect, 
we did not expect to find a significantly different value. Indeed 
we obtained: Mp cor = 5.807338 x 10°Mo, consistent with the 
previous one. From the virial black hole mass and by using the 
same value of bolometric luminosity reported in Sect. we 
computed an Eddington ratio of 0.86*}35, which is consistent 
with that computed from the Ojjn¢. 


Mgu = 10°°° x ( x ( yr (AD) 
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